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XNTROOUCnON 

It ha« long been recognized that the top few ffllllloeters 
of a surface exposed to space represents a physically and 
chenically active zone with properties mich different from 
those of a surface in the environo^nt of a planetary ataosphere* 
It is in this zone that the naterial comprising the surface of 
the moon. Mercury, the asteroids, and similar bodies achieves 
physical and chemical maturity as a result of exposure to the 
solar wind (SW), solar flare (SF) particles, galactic cosmic 
rays (CCR), heating by solar radiation, and to the bombardment 
by meteoroids of a vast size range. 

The recent work of several groups continues to clarify the 
nature of the processes involved at space-exposed surfaces. 

These include ion implantation and redistribution (Tamhane 
et al „ , 1979; Warhaut ^ 1979) and sputtering effects (see 

review paper by Pillinger, 1979), the origin and distribution of 
metal particles determining the ferromagnetic resonance (FMR) 
maturity index (Morris, 1980), shock effects and the formation 
of giacL and agglutinates (Schaal and liBrz, 1980), and tlie 
effect of intergrain adhesive forces (Housley, 1980). 

There has been a clear need for a quantitative synthesis of 
the various processes into a self-consistent picture of surface 
evolution. This synthesis is especially Important as different 
experimental approaches have often resulted in apparently 
conflicting interpretations of soil history ;it specific sites. 
Recently Wleler (1980) have attacked the problem experi- 

mentally by measuring noble gas concentrations in Individual 
grains and comparing these with track distributions from the 
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same soil sample. This approach la vary aaaful as it is often 
risky to correlate data neaaurad on diffarabt s a apl a a . 

The present work represents a step in our continuing 
effort to achieve a model synthesis. This study draws upon 
the experience gained from earlier models by the author and 
other groups ( e^>£,. » Gault at al. , 1974} Arnold. 1975; Duraud 
gt aJL., 1975; BorggLAl*. 1976; Comstock, 1977, 1978). 

To obtain the results described here developed a new 
set of computer programs called MESS.2 (Comstock, 1979, 1981). 
Ihe design of the MESS.2 system Is described in our Final Sep<MTt 
for NASA Contract NAS 9-15582 submitted September 25, 1979. 

This current model differs from earlier work in that surface 
processes are broken down as a function of size scale and 
treated in three distensions with good resolution on each scale. 
Since the smaller size scales are on the order of a grain 
dlauieter this approach derives explicitly the tiiiM and depth 
dependence of parameters &s functions of grain size as well as 
of sample size, mixing depth, and surface layer residence time. 
Factors relevant to grain size dependence have been discussed by 
several authors (Eugster aJL* » 1975; Criswell, 1975; Becker, 
1977; Bogard, 1977; Criswell and Basu, 1978; Housley, 1980). 

Our results here apply to the development of soil residing 
near the surface for the first time and hence should describe 
immature to submature soils. The development of mature soils 
through further exposure and mixing can also be derived using 
these results. The current model is based on lunar conditions 
where considerable data exist to provide guidance and calibra- 
tion. The model parameters can readily be adjusted to describe 
asteroid regollths und other space-exposed bodies. 
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mss *2 HO&EL 

Grain redistribution* 

MESS *2 breaks the surface soil dovn to a series of 
relevant size scales or reglsMss, siaulates each regioe separate** 
ly and combines the results in an appropriate statistical 
manner* The structure of a given regime is shown schematically 
in figure 1* The depth scale is determined by the bottom of the 
model layer » depth D^* A surface zone of depth is defined 
corresponding to the SK>del layer of the next smaller size 
regime* Similarly D corresponds to the surface zone of the 
next larger size regiiM!* 

For a given regime a range of Impact event sizes is 
simulated such that the largest events will generally just 
disrupt the full model layer depth, while the smallest events 
will generally just disrupt the surface zone. Hence the events 
occuring in larger size regimes will determine the local resi- 
dence time for the model layer while events considered in the 
given regime control the local residence tlioe of the surface 
zone. Events in the given regime also primarily control the 
ejection of material into the surface zone and subsequent 
reburial below it over a time scale comparable to the average 
residence time of the model layer. Details of the soil history 
within the surface zone are controlled by smaller size regimes* 
The events that effect a given regime will create a mixed zone 
of variable depth D which containr all material ejected by the 

K 

simulated events* 

Clearly this procedure only approximates the true behavior 
of the soil, we BXist emphasize that our purpose Is not to 
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attenpt to slimlate the soil precisely but tmly to efficiently 
nodel the physically Important processes goveiaing Its observed 
characteristics. 

In the current work we use regimes. 

In each regime the smallest crater slimilated has a diameter 
equal to most events have a diameter lesa than about 3 0^^. 
Depths are monitored on a continumis scale while the surface is 
divided Into an heM;p‘ .tal grid, to enhance inherent circular 
symetries, of about 3(KK) cells where the cell diameter is no 

■/ 

larger than D^. This assures good resolution for all events. 

The location of up to 12,000 iivlividual particles of a given 
grain size is monitored by 1ESS.2. For about 200 of these we 
also BK>nitor central track density and surface-zone residence 
time, initially these tagged grains are pristine and randomly 
distributed throughout the model layer. 

For each grain size there will be a size regime such that 
the surface zone corresponds to a grain diameter and the smal- 
lest event will barely eject some grains intact. This is the 
smallest size regime modelled fur a representative grain size, 
referred to as the base regime. In this reginte the individual 
particles move in a matrix of order-of-imjgnetude finer material 
so as to approximate a typical size distribution. In this way 
the graininess of the soil and shielding of individual grains 
near the surface is taken info account in deriving distributions 
of exposure tinx* and track density. A detailed size distribu- 
tion of the matrix is not considered, so possible shielding by a 
thin layer of relatively very small particles is included only 
thriaigh a constant exposure f«vt«>r (f) whlcti Js discussed in 
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detail latar. In larger else regiMS all grain aliM (ei^ to 
behave alike * aaaller grains iMve saallcr else base reglMis 

and therefore generally experience a dlffei^nt dyiiMSlcal history 
near the surface. This contributes to the overall grain else 
dependence of characteristics requiring surface exposure. 

We are interested In grain sizes ranging from less than 10 
aicrons to greater than 1 nua In dlaneter. Soil saoples usually 
average over intervals ranging froa 1 ma to 1 ca In soil colums 
to several centiaeters in scoop saaples. Therefore In this irork 
we investigate the size regimes “ 100 microns (base region 
for lO micron grains), ■ 1 an (base reglsM for 100 micron 
grains), ■ 1 cm (base regime for 1 wbl particles), and ■ 

10 cm. The results will apply to soil which has resided once 
within 10 cm of the surface and should characterize immature to 
suimiature soils. 

The motion of ejected uaterial is followed in three dimen** 
sions during each Impact event according to an ejecta origin 
model. The ejecta map we have developed for MESS. 2 is shown 
in figure 2. The map was designed to be consistent with the 
results of StSffler e^ (1^73), one of the few experimental 
studies which clearly defines the ejecta origin for soil craters 
in the 10 cm diameter range. The work of Vedder (1972) suggests 
that the ejecta distribution for craters on the order of 100 
microns diameter is probably much broader than that shown in 
figure 2. This will not be important to our results htwever 
since on this small scale our individual grains are deposited in 
essentially a single layer after ejection and hence their exact 
location is not statistically significant. For craters larger 
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than a few centiaetera ejecta sapping Md atratlgraphy reveraal 
occur In HESS • 2 eventa* 

Our auboeter cratera are baaed on the few experiaental 

Btudlea available Breslau, 1970; Vedder, 1972; StHffler 

et al. , 1975; Gault and Wedekind, 1977). The work of Gault 

and Wedekind (1977) can be applied to the moon becauae they 

determine the dependence of crater paraaetere on gravitational 

acceleration. They find a crater depth H to rla dlaaeter 

ratio H/D- ■ 0.2 for a 15 ca crater, easetitlally Independent of 

gravity. Stuffier et al. (1975) obtain similar results but show 

further that the lower part of the crater is not ejected but 

compressed Into the crater floor. Their craters In weak quartz 

sand show an ejected depth to diameter ratio closer to 

"0.1; In more cohesive lunar soil the ratio may somewhat 

greater than this. Vedder (1972) investigited much smaller 

craters of about 100 micron diameter in a nore cohesive test 

soil and obtained H/U " 0.2 <- 0.3 much of which represented 

compressed material. We have adopted a va.uc of iU/D„ " 0.15 

over the full size range studied and only laterial actually 

ejected is considered redistributed. MESS. 2 ignores soil 

compression and density variations in general as second**urder 

effects on the irradiation history of a grain. Depth scales 

3 

In this paper apply for a soil density of 1.5 - 2.0 g/cm . 

We have taken local topography Into account such that crater 
volume is maintained while generating a smooth crater that 
follows local slopes. 

The sizf.‘~f requeue y distribution of small soil craters 
cannot be determined bv dire* t count as the lunar surface is 


QUAUTY 


8 


dpmu^iMjnr 


•ttturaced with craters on this ecele* Xneteed we itee the else* 
frequency dietrlbuti<Hi of pite oheerved in rock eurfecea end 
tranaletc thia to aoil cratcra using the pit slse to projectile 
MSS relation and aastuMd velocity givm by SVrs et,^* (1975)* 
For pit diaseter DP in ca and projectile •imtgf E In erga ttelr 
calibration yields Dp " 2«lxl0~^ The work of Gault a^ 

Wedekiiwl (1977) would then indicate that a lunar soil crater of 
D^ ■ 25 CB would correspond to a 1 ca rock pit, or * 25* 

However Gault and Wedekind used dry quarts sand to alnlsdse 
strength effects so we should assuoe that 25 represents a 
aaxiaua value for this ratio In the lunar case* Braalau (1970) 
obtained results slailar to those of Gault and U^eklnd. Vedder 
(1972) finds that soil craters are saaller in mre cohesive soil 
Indicating that strength effects are l^ortant in this slse 
range* Gault and Wedekind reach a slailar conclusion even in 
the case of dry quarts sand* The results of Vedder (1972), 
using a much lower projectile energy and a more cohesive soil 
indicate that a similar ratio D^/Dp * 10 - 20 should be valid 
down to 100-micron diameter craters* 

For the size-frequency distribution of Impact pits we can 
CMibine the results of Httrz ejt al* (1974, 1975) aiul Morrison and 
Zinner (1977) to obtain an approxiaate analytic fora: 

(1) NOD^) - A ® 

2 

where N is the cumilative pit production rate per cm -My, DC is 
/Dp, and the constant 0.005 reflects a jog in the 
pit distribi'.tlon at 10 - 100 microns. A is about B/1000 but the 
first term in equation 1 is not important for the redlstrik>ution 
of grains with diameters of about 10 microns or more* 
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Il8rz ^ £jL. (197A) use B - 8.3x10 Comstock (1978) suggested 

that a value up to 10 times higher might be more consistent with 

track gradients in dusty rocks. HBrz ££ £l. (1975) discuss 

sources of uncertainty in the long-term micrometeoroid flux. 

Considering the uncertainties in D /D and in meteoroid 

flux we have defined a model time denoted MT such that I MT 

corresponds to 1 Hy (10^ years) in the case of 0^/Dp - 20 and 
-A 

B ■ 8.3x10 . The value of MT is Inversely proportional to 

production rate N and, for D i 1 cm, Inversely proportional to 

3 

(DC/DP) , so the uncertainty in the time scale MT might be as 
much as an order of magnet'ide. However we show later that an 
assignment of 1 MT ■ 1 My is supported by calibrating the model 
results against lunar data. Note that the time scale depends on 
the cohesiveness of the soil through the ratio Dj,/Dp. 

Grain alteration. 

In addition to the redistribution of soil particles as 
discussed above particles also experience the effects of 
irradiation, erosion, f lagr.ientat Ion, aggregation, and sliock 
metaiiiorphosis . 

MESS. 2 accumulates particle tracks continuously in the 

tagged grains using the solar flare track production rate given 

hv Comstock (1978) of 7x10^ d cm ^-My ^ for depth d in 
2 

g/cm below a plane surface, and the galactic cosmic ray track 

rate given by Walker and Yuhas (1971) which predominates below 
2 

about d ■ 0.5 g/cra (about 3 mn of soil). In the base regime 
for a given grain size tlie central track density is acciunulated 
at every depth the grain resl«les. In larger size regimes tracks 
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are accuimilatedi only at daptha balov ^ Mirfa^ mmim aiM i^^d 
to those accuBuIated in Mailer regiMa aa a funetlM of aarfaM 
zone realdence tine* In thia way tiie CMtrlbitti^ frM Mch 
level of raaidcnce la obtained ei^licitly* 

Tte neaaured flux of aolar wind l<ma la revletied by 
Pillinger (1979). The obaerved concontration of implanted SW 
Iona dependa on retention charactcrlatica different for each 
apeciea and boat mineral aa well aa aurface expoaure time* In 
this paper we derive only with eiqpoaure tine and raaldMce tlmas 
within each depth zone* In tlM base regime tlESS*2 monitors 
the actual surface exposure time of tagged grains* In larger 
regitsea the time spent within the aurface amie is mmltored* 

We are concerned here with processes which effect the 
irradiation record of an enserable of grains over a time scale 
appropriate to a single residence episode within 10 cm of the 
surface* We are particularly concerned with effects important 
at the surface or varying with grain size, ^K>ng these are 
direct erosion of exposed grains and fragmentation of grains by 
direct impact. We find that these strongly effect the oagMtude 
and grain size dependence of the irradiation record* 

We assume that a grain of diameter will be fragmented 
if it receives a direct impact generating a pit of diameter 
On > D^/8. This value is uncertain and varies with condition 

r U 

of the grain, but smaller pits have been observed on soil grains 

vdiile a pit approaching must certainly fragment it* For the 

o 

-4 

microsMteoroid flux used to obtain equation 1 with B • 8*3x10 
this yields a lifetime against fragmentation for 10 micron, 
100 micron, and 1 nan grains of about 3 My, 2 My, and 1 My, re** 
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spectlvely, with conilderable uMtrtaintf. With 10 tlMs higher 
flux these becoac 0.3 ffy, 0,2 tty, end 0.1 tty, respectleely. We 
have perforaed our cslculatltms for lK>th of tlMse sets of velties 
but since the shorter llfetlaes hsim auch M>re effect t for 
clarity we coaparc only the results for the shorter gainst 
the results assuming no fra^ientatlon. 

FragiMnted grains are re»>ved froa their site population. 
Surface fragmentation will tend not to effect the low track 
density paraaeters such as the mlnlaiM and quartlle since It 
will reaove only those grains which have resided on the surface 
long enough to attain the highest densities. It can have a 
strong effect on solar wind concentration however since It 
removes those particles which have received the greatest SW 
fluence and hence control the total concentration. 

The erosion of exposed grain surfaces by mlcrce^teorolds 
too small to frajpaent the grain is generally not a significant 
process because of the low flux Involved. The most important 
source of grain erosion relevant to SW ion implantation is 
atomic sputtering by solar ions (Fillinger, 1979; Flavlll et 
j^.yl9B0) who report a consensus of recent measurements suggest- 
ing a rate of 0.U2S to U.045 2/year at the lunar equator. For 
a typical Implantation depth ot 3(K) 2 (Fillinger. 1979) this 
leads to an Implantation/erosion equilibrium on a time scale of 
T^ “ 0.01 Hy. We present model calculations both with and 
without an effective SW saturation with this time constant. 

This saturation time depends on implantation depth and so 
will vary with species (Plliliij’.er . 1979). It has been suggested 
by Signer . (1977) that diffuslonai processes also leads 


to a certain saturation tine constant. ;^ialar at al. <1980) 
note that diffusion and sputtar eontrollad saturation uill have 
siailar effects. Likewise the MESS.2 calculatioos depend onljr 
the tine constant, whether this tine constant arises fron sput- 
tering or dlffusional processes is tMt readily deterninahle 
using our results. 

The fomation of agglutinates by shock conpaction during 
nlcroneteorotd iapact also occurs within the top ailllasters. 
This process loads the agglutinate with irradiated grains and 
reduction !oetal. Tims residence tine within the first few 
mllliiaeters controls the noble gas and FMR characteristics of 
agglutinates. MESS.2 is designed to n<mitor glass and aggluti- 
nate production but these id.ll not be discussed in detail hare. 
Class and agglutinate production and shock cosstinution through- 
out the soil volune becooe increasingly inportant as the soil 
natures and these processes mist be taken into account 
explicitly in any simulations of nature soil. 
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HODEL RESULTS 

Hont«~C«rlo calculatiofit. 

To illuatrato th« roaulta produced by the MESS*^ aodol tm 
•how in figures 3-14 sMie of the residence tl«e end track 
density diatributione derived ^'or lOCHaicron dleaeter greina. 
Soil greina with dleaeter on tlie order of 100 alcrona ere typi- 
cally the more coaron aiaed pe ‘tides in lunar aolla* Siailer 
diatri^tiona were elao derived for repnaentetive particle 
dimetera D. ■ 10 aicrona end 1 mb. Figuzea 3*9 indicate 

V 

diatributione not yet including f re^enteti<Hi loaa. In figures 

3 * 8, end 9 the date points refer to a typical coaputer run. 

The base regime for 100*micron grains has a aodel layt*r of 

thickness 1 mb, identified as the regiae. All craters 

of dlaoMter 1 bm < < 1 cm are considered, the saaller ones 

greatly predcHslnating. Figure 3 shews the derived correlation 

between central track density and surface-exposure tine after a 

typical layer residence time " 0.1 MT. Recall that model 

time MT aost probably is equivalent to 1 My. The rectangle on 

the left represents the range of track densities for grains not 

yet exposed at the surface. The solid diagonal line represents 

the track density a grain would Itave if it accuaulated tracks 

only while it was on the surface for a time t, while the dashed 

line represents the maximun density possible after time 

thus enclosing the allowed zone. All grains with track density 

on the order of 10 cm or mure have resided directly on the 

surfa< All exposed grains have a central track density of at 
<2 

leat cm and will tiave strong track density gradients. 

During a typical regls^ residence time of about (>.) MT 
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nearly half of the grains will have hMO ei^aed at the surface. 
For these » cuoulative ex{K>Sttre tiae diatrihutlona are ahoiin la 
figure 4 for different resid^ice times T^. Although only half 
of the grains in the 1 bmi layer have been exposed the average 
expoETure tiiM is about 0.1 T , indicated the effect of shield* 

nn 

ing by fii^r material. Note that shielding by larger grains <m 
the order of 1 bob or more is excluded by the constraint that the 
simulated grains be within 1 mm of the surface, but will be 
accounted for in the next larger regime. 

Figure 5 shows some track density parameters resulting from 
residence within the top 1 ram. Slopes greater than 1 in figure 
5 reflect variation in a grain's depth due to Impact gardening 
together with the steepness of the SF track production profile. 
Since grains which have a measurable track gradient as defined 
by Comstock (1977) muse have resided within about 1 mm the para- 
meters in figure 5 refer to that component of the central 
track density which would be related to an observed gradient. 

The motion of 100-raicron diameter grains is also followed 
in the next two larger regimes, for • 1 cm denoted the CM 
regime and for D. • 10 cm denoted the DM regime. For the CM 
regime we show in figure 6 some representative distributions of 
non-sero residence time within the 1 mm surface zone, identified 
with for two typical 1 cm layer residence times Only 

about a third of all grains have resided within 1 nun and these 
for an average time of about 

In figure 7 we plot parameters describing the tracks 
accumulated during residence between 1 mm and 1 cm depths, this 
represents the major contribution to central track density which 
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does not result in s track gradient in the grain* The Dlnima 
density refers mainly to grains that have remained near 1 cm* 

For those grains which have been cycled through the upper 1 mm 
surface sone the track densities Indicated in figure 7 represent 
only part of the total accumulation and naist be added to the 
tracks acquired within 1 mm (figure 5). Mostly these are grains 
which have resided in the upper part of the 1 cm layer and hence 
have higher partial densities the median for this sub- 

set of well-irradiated grains is given by the dashed line in 
figure 7* Note that the maximum contribution to the total track 
density from residence below 1 mm (figure 7) tends to be compar- 
able to the first quartlle density contributed during residence 
within 1 mm (figure 5). 

For the 10 cm (DM) regime figure 8 shows typical distribu- 
tions for residence time T > 0 within the top 1 cm for two 

CM 

residence times T ^ within 10 era. The average values are about 
DM 

Tjjj^j/3. The partial track density accumulated below 1 cm is 
obtained and taken into account, however it represents a factor 
of about ten less contribution to the total density than does 
residence between 1 mm and 1 cm. 

Since most lunar soil sample intervals are centimeter or 
subcentimeter in size we are concerned more with the structure 
within the 10 cm regime, distinguishing samples in the mixed 
zone from those below the mixed zone. Some properties ot the 
mixed zone are shown in figure 9 for the DM regime, similar 
properties hold in all regimes. In the botton part of figure 9 
is plotted the cucmlatlve distribution of depths disturbed by 
impact events after a 10 cm layer residence time of 2 MT. In 
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all our plots the depth la assured raletive to the cia^rent, 
olMervable surface not the initial surface* so It Includes sites 
tdilch have been recently excavated below the previously die* 
turbed level and have only a thin mixed sone deposited siiu:e« 
This distribution was collected over a surface area of radium 
about 25 cm* a single core site would not reflect the broad 
spread even though it occurs quite locally* 

In the upper part of figure 9 we have plotted the correla-* 
tion between depth and residence time within the surface seme 
for Uie srae layer residence time and surface area as in the 
lower part. This plot makes clear that soil properties which 
depend on surface zone residence will be distributed equally 
throughout the mixed sone and that no general c>' .relation with 
depth is to be expected. 

This does not mean that the mixed zone will be found to be 

uniform at a particular core site* because the data in figure 9 

were collected over an area much wider than the area of a 

sampled soil column. We find that over a typical core sample 

area there can be strong variations with depth. The variability 

of centimeter-interval samples In a soil column will reflect the 

variability of the Q1 regime as a function of T_ . allowed by 

the broad distribution shown in figure 8. Moveover there may be 

considerable variability on the millimeter scale reflecting the 

variability of the til regime as a function of T * also broadly 

rin 

distributed (figure 6). Evidence for strong variation of track 
densities over a few tm has been reported (£>£* * Goswaml et * 
1976). 

MESS. 2 results show that ram- and cia-scale variability can 
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«&8ily b« very local in origin* WImo eueh variation in the 
exposure and irradiation record iv observed experiiMntally it 
shcxild not in itself be considered sufficient evidence for a 
time variation in charged particle or nicroaeteoroid flux or 
for a different soil origin unless a definite correlation with 
depth can be established* 

The time develointent of the median depth of the mixed seme 
is sho«m in figure 10 for each regime by the solid curves, which 
were collected over a area of radius 2*5 D^* A regime was 
defined as a layer of depth D, which was assumed to avoid local 
obliteration by large impact events for a time T^* The median 
depth mixed by all events over a sufficient!/ wide area will be 
Just the envelope of these curves indicated by the dashed line* 
The more one of the solid curves deviates from the dashed curve, 
the less likely it becomes that such a shallow depth will occur 
except very locally. Two things must be kept in mind when using 
plots such as these. First the mixed depth is broadly distrl" 
buted, as indicated in the bottom part of figure 9, and may 
c<xnmonly be at least a factor of 2 ..a either direction from the 
median shown in figure 10. Secondly over a very small area such 
as is sampled by a soil core the mixed depth can vary erratical- 
ly with time (either up or down) and will not necessarily follow 
in detail any of the curves shown, except statistically. 

In an earlier model (Comstock, 1977, 1978) we assumed that 
all grains had resided within 1 mm of the surface and followed 
their motion as they were buried through the top few millimeters 
over a time of about 1 My which was called a surface exposure 
episode (SEE). The scenario assumed then is therefore equiva- 
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lent to grains residing in the nixed sone of <Nir CM regiae aodel 
which we see frost figure 9 can be mixed to a few milliiMters by 
local events in a time 1 MX* Although the current aodel treats 
detailed motion of the grains more accurately our earlier work 
is still approximately valid as a subset of our current results 
and probably represents a coamon situation. The current model 
covers a much broader range of soil residence conditions as well 
as deriving exposure and residence time distributions. 


Sample synthesis-tracks. 

The results obtained for each individual regime must now be 
combined to obtain the characteristics of a realistic soil 
sample containing grains which may have resided in more than one 
regime. We discuss first the distribution of particle track 
densities for lOO'-mlcron diameter grains which is representative 
of the grain size us’ lly used for track studies. 

The partial track densities accumulated during each regime 
as a function of residence time are first added together accord" 
Ing to the derived distributions of residence times. This 
results in a series of time-dependent relations for the track 
parameters as a function of the shallowest regime in which the 
entire sample has resided. These relations for the minlnxim, 
first quartlle, and median track densities are shown in figure 
11. The time scales are lined up so average values correspond. 
It should be noted that since all the track distributions in 
this study are derived using a finite number of grains the mini- 
mum density in fact corresponds to about the first percentile. 

A similar caveat should be made for all experimental minimum 
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track danaitiea. The true theoretical BlnimuB density will 
always be soDeidiat leas and not likely to be observed. 

In figure 11 the curves i&arked Ifi describe a sample which 
has resided in the top millimeter but includes the contribution 
from residence below 1 nmi. The curves should also apply well 
to any no-sized sample taken from the mixed zone of a 1 cm surf- 
ace layer. For mm-samples taken below this there tsay be an ad- 
mixture of less-irradiated soil. Similarly the curves marked Qi 
apply to cn-sized samples taken from a mixed zone up to several 
centimeters thick and the DM curves apply to dm-sized samples 
taken from a mixed zone several decimeters thick. 

The data which lead to figure 11 can be presented another 
way in terms of a time vs. depth plot for a particular parameter 
as shown in figure 12 with the first quartile track density 
indicated by contour lines. The dashed line from figure 10 
giving the overall mean mixed depth is also indicated by the 
daslied line in figure 12. The contour lines in the unmixed 
region reflect the shape of the track production profile in 
spite of some depth smearing due to gardening of the soil above. 
Witliln the mixed zone the contour lines become distorted and for 
mixed depths on t!ie order of centimeters the very shallow depth 
dependence of the average quartile density is apparent. 

The curves In figure 12 were derived assuming a statistic- 
ally average sample, at a given sample site the contour lines 
within the mixed zone will depend on the actual series of 
impacts effecting the site sampled and on local mixed depth. 

Me can see from figure 11 that any depth dependence that might 
have been expected within the mixed zone can be readily masked 
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by the wide variation in T_ and alloiied throi^hout the 

nil Un 

mixed zone, refer for example to the upper part of figure 9. 

Only within the uppermost 1 m of soil is a depth dependence 

likely to become apparent because of the steep track production 

profile there* An Interesting feature in figure 12 is the 

suggestion that the quartlle track discontinuity at the bott<MS 

of the mixed zone becomes sharper as the zone becomes deeper* 

In order to help calibrate our model time scale MT we have 

Included in figure 12 a data point. Indicated by the circled 

dot, obtained from the top of double drive tube 74002 collected 

at Shorty Crater* This is one of the clearest examples in the 

literature of an apparently fresh soil residing within the top 

10 cm for the first time and, according to Crozaz (1979), is 

mixed to a depth of 6 cm over a time of 10 My. Crozaz finds a 

6 

quartlle track density within the mixed zone of about 3x10 
-2 

cm . Assuming that the Shorty Crater data represents an aver- 
age sample it is consisent to assign 1 MT ■ 1 My, thus agreeing 
with the quartlle track density shown in figure 12 as well as 
with our estimate above based on the crater production rate. 
This time scale also allows our mixed depths (dashed lines in 

figures 10 and 12) to match smoothly with mixed depths deduced 
22 2 ^ 

from Na and A1 profiles in lunar cores (Fruchter et al. , 

1976, 1978) and from I /FeO profiles (Morris, 1978)* 

s 
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Saaple synthesii-»‘«xpoaure cIim. 

( We turo now to the derivation of average reaidence-'tlve 

functlona which would allow the calculation of other irradiation 
effects, aolar wind iaplantation, FMR, etc. We are intereated 
in the concentration per unit maaa of aooe given species 1 

introduced at a rate per unit area per unit tiaie into soil 

grains residing either at the surface, denoted EX, or elsewhere 
within reglne Ol, etc . We wish to obtain as a function 
of residence tine T, within a soil layer of depth D, . The 
desired concentration can be written in the following form: 

f"* 

(2) Cj(T^) - Fj I (2f/pDg) <tag/M. 

In this relation dm./M is the differential mass fraction for 

G 

grains of diameter D.. The factor 2f/pD_ gives the grain area 
_ G G 

( / exposed per imlt mass for grains of diameter where p is the 

grain density and f is the exposed area of a grain divided by 
its cross sectional area. The function represents the 

effective average time spent within any zone X as a function 
of layer residence time Tj^ for grains of diameter D^, where X 
and L nay stand for EX, Ql, DM, etc , and the subscript G 
stands for in centimeters. We use represent 

grains knoim to have come from the nixed zone of the layer L. 

We are concerned here with the behavior of the model inde- 
pendently of the rates and the exact distribution of grain 
sizes. Therefore we calculate effective area-time production 

functions C , for each grain size D , defined as follows: 

X o 

(3a) ■ (2f/p)[Ijj(T^)g/DgJ for the whole layer L and 

(3b) G*(Tj^)g “ (2f/p)II*(Tj^)j,/Dj,J for the mixed zone of L. 

We assuiM that f Ideally is on the order of 1 and that 2f/p is 
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on the order of 1 in cgs units. Therefore in figures 13 to 18 

the appropriate grain diaaeter in cn. In fact p will vary with 
grain composition and f will depend on grain shape, orientation, 
and partial shielding. Enhanced partial shielding by very saall 
grains sticking to larger ones nay be particularly Important 
(^.riswell and Basu, 1978; Housley, 1980) so it must be kept in 
mind that much uncertainty and some additional grain size depen- 
dence may still be burled in the factor f, as we discuss later. 

Before we calculate the functions we need to specify 
for each regime the relations for the distribution of non-zero 
surface-zone residence times and for the fraction of grains 
which have resided in the surface zone as a function of regime 
residence time. We will also need mechanisms for modifying the 
exposure time distributions to take into account the effective 
surface saturation of an accumulating species and the loss of 
grains by impact fragmentation at the surface. 

From the MESS. 2 Monte-Carlo simulations (je , figures 4, 

6, 8) we find that in each regime the differential distribution 
of non-zero surface-zone residence times is described well by: 

(4) - [exp(-T2/K^Tj^)J/(K^Tj^)(l-exp(-Kj/K^)l 

for 0 < Tj < KjTj^ .na - 0 for > KjT^. 

where is the main parameter, is a cutoff usually about 1. 
In a base regime T stands for the actual space-exposure time 
t. The fraction of grains which have been cycled through the 
surface zone is approximately given for typical T by the form: 

(5) f(T)_-f (l-exp(-T./T^)J for T^/3 i T < 3T,. 

eLGm Lr r Lr 

For the fraction f^ will continue to increase slowly 


I* can be recovered from and Gj^ by multiplying by 
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but this h«8 little effect on the reeulte* Pereaietere f end 

ti 

define a nixing efficiency and tine c<metent| reepectlvely. 

The values found for the various constants are given In table 1 
for each regime. Note that grains residing In their base regime 
behave somewhat differently than do smaller grains In that 
regime, due tu the gralnlness of the model soil* 

The statistical method outlined In the following paragraphs 
represents an approximation to the actual grain depth history 
within a surface cone, MESS* 2 was originally designed to follow 
more precisely the grain depths In each surface zone by using 
smaller scale statistics continuously during the computer run 
for a given reglTO* This was found to be very time consuming 
without giving significantly different results than the methods 
described in this paper. The effect was to alter slightly the 
detailed shapes of the exposure and track distributions but not 
the average exposure tine or the maximum and minimum Irradiation 
levels; the additional complexity is not justified In view of 
current uncertainties in other model parameters. 

We now derive the average residence time functions 1^^ in 
the following manner. For a given grain size we first start 
with the base regime, for lUO-nicron diameter grains this is the 
IIM regime. Since the desired characteristics of grains anywhere 
within the nixed zone reflect residence in the surface zone we 
first average over possible surface-zone residences. The aver- 
age space-exposure time for grains from the mixed zone of the FM 
regime in the absence of fragmentation and saturation is: 

r 


If we sample the entire I mia lay«>r then we must add the unmixed 
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grains as well by renormalizing the Integral aa follows: 

In equations 6 and 7 the functions f and f use paraneter 

8 « 

values from table 1 for lOO-mlcron grains In the HM reglM. 

Now we must correct this to take Into accoimt fragsientatlon 
and saturation effects. We assume that the fraction of grains 
surviving loss by fragmentation for a surface time t Is: 

(8) f^(t)jj - exp(-t/Tj(D^j)). 

We approximate grain-surface saturation processes by multiplying 
the integrand in equation 6 by a saturation function given by: 


(9) fj.(t) - (Tg/t)[l-exp(-t/T^)]. 

The functions and I* will then represent effective average 
times with saturation built in unless goes to infinity. 

With the correction factors for grain fragmentation and 
possible saturation included the average exporure tln» in 
equation 6 for the mixed zone of the top 1 mm bec(xaes: 


( 10 ) 




ami 


HH 








where F 


r^MM 


Equation 7 for the whole 1 mm layer becomes: 

where f i? f (T^) and the fragmented grains have been 
e e nn .ui 

removed from the 100-mlcron grain size population. 

We now construct the effective average exposure time for 
lOO-mlcron grains that reside in the upper 1 cm of soil, 
simulated by the CM regime. Fur grains from the mixed zone of 
a 1 cm lajer we ideally would average over 1 mm surface-zone 
residence episodes as follows (suppressing subscripts EX and .01 


for clarity; ) 
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(12.) I*(T^) - p, 




y:i 








where the first tern is contributed by that fraction of grains 

which have resided only once within the 1 nm surface zone, 

Che second term is contributed by that fraction of grains 

)4iich have resided twice and only twice within the surface zone, 

and so on for high.i;r order terms. The distribution f' gives 

6 

the frequency of surface-zone residence times for individual 
episodes within the surface zone, as a function of layer time. 

In two independent cases equation 12a will reduce to the 


much simpler form given by: 




- f 

J 0 




where f is the total surface-zonu residence lime distribution 
g 

such as that shown in figure 6. The two cases occur under the 
following conditions, either: 

(a) p "0 for all n > I, so that p ■ 1 and f' - f ; or else 
n 1 , g g 

Is proportional to T„,,, so I(T' )+I(T-T' )-l(T) 
fcX MM liM 

which along with the definition of f in terms of f' gives 

g g 

equation 12b from 12a* 

Our MESS.2 Honte-Carlo calculations give the distributions 

f' and fractions p so equations 1 2a and l^b can be ct)mpared. 
g *n ^ 

We find that ft^r low values of the condition foi case (a) 
are generally satisfied (but not case h). Fur higher case 
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(b) !• ntarly Mtiafiad in th« abtenca of grain fragponcation 
or aacuration, although two or thraa aurfaco-toiM aplaodaa 
bacMM coanon; on tha othar hand if bacoMB aaturatad 

than aultipla aurfaca-iona apiaodaa bacoaa inaffactiva again dtM 
to loaaaa or grain*-aurfaca aaturation* and tha co^ltlon for 
caaa (a) la affactivaly approached inataad* In all eaaaa 
aquation 12b baconaa a uaaful approxiaation, 

VlMn tha entire 1 ca layer is aaapled ua will have: 


(13) “ V^CM^Ol 

If we need tha a^rage tiae spent W4.iUin 1 aa as a function of 
Trnjt aaaplas froa the alxad s<ma of tha 1 ca layer, we use: 

laM 




.r 

.01 J 


CM 


^MM ^g^^lM’^CM^Ol **^MM' 


For the whole layer we use: 


^^5) ^mm^^cm\oi 


^e^^CM^Ol ^ii^^CM^.Or 


In equations 12b through IS tha functimis f and f use tha 

8 • 

parameters in table 1 for lOO^icron grains in the Q1 regime. 

Sinllar relations can be obtained for saaplea from a 10 cm 
layer. The most useful of these give average space-exposure 
time and average Ina-layer residence time for samples froa the 
alxad sons of the 10 cm layer approximated respectively by: 


(16) 


^EX^^CM^Ol ^g^^CM’W.Ol 


(17) 


^W.Ol " J Q 
C ^DM 

^MM^^DM^.Ol "Jo 


01 


where the parameters froa table 1 for lOO-tslcron grains In the 
DM regiiM are used In f^. Equations 16 and 17 are used for the 
same reasons discussed in reference to equations 12a and 12b. 
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If th« functlMia found abovn ar« now dluidnd by grain 
dlaaeter • 0*01 ca w« obtain tha aff active araa-tlaa piroduc~ 
tlon functlona aa diacuaaad aarliar. Production functlona 
wvira derived aiailarly for reprcaentatlve grain diaaatara ■ 
0.001 and 0*1 ca* boae of these reaulta are shown in figures 13 
tj 18* In all of these figures the curves labelled A, B, or C 
rjfrr to different aasuaptlons for the fragVMntaticHi tiaa const- 
ant T, and the saturation tiae constant T aa defined earlier* 
Tite curves labelled A were derived aasuaing no fragaentation or 
s-ituratlon. The curves labelled B assiae only grain fragmenta- 
tion, using values for of 0.3, 0*2, 0.1 MT for 10-aicron, 
l(K)-aicron, and 1 wa diameter gmiris respectively (except in 
figure 18 discussed later). With an identification of 1 MT ■■ 

1 My these correspond shorter fragaentation times dlscus- 

8< d earlier, 'ive curves libelled C assume both these saoie frag- 
mt ntation losses; plus n surface saturation time of T^ • 0.01 MT 
(except in flv' re 5 8) again corresponding to a plausible value 
discussed ibo^e for 1 MT ■ 1 My. 

?lgures 13 to 15 give the effective production functions 
fur direct space exposure. In figure 13 we plot 
d«*scribing aa-slzed samples taken from the mixed zone of a 1 cm 
surface layer, figure 14 shows describing the entire 

1 cm layer including uiiexposed grains below the mixed zone, and 
figure 15 shows describing a 1 cm sample from the 

-<xed zone of a 10 cm layer, all as functions of the whole layer 
reaidence time. Figure lb gives ) , the production 

function for residence within I rim fur gr.iins In a 1 cm sample 


t.iken from the nixetl zuiu* of a lu cm layer 
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DISCUSSION 

Figures 13-10 all have Important features In connon* In 
the absence of fragmentation or saturation (curves A) 100-mlcron 
dlao^ter grains show a lower production function than either 10- 
mlcron or 1 mm diameter grains. This results from a combination 
of effects: A larger fraction of grains (t'ur model matrix) Is 

available to lightly cover 100-mlcron grains compared to 10- 
mlcron grains In their respective base regimes, ^.js. very close 
to the surface. 1 mm grains would be even more readily covered 
by finer material but the steeper size-frequency distribution of 
craters ( eg rre spending to >500-mlcron diameter rock pits) in the 
CM base regime for 1 mm grains gives them a better chance to be 
exposed by shallow events while residing very near the surface. 

In the presence of grain fragmentation (curves B) the 
effect of fragmentation by direct impact while on the surface is 
much stronger for 1 mm grains due to longer surface exposure but 
similar fragmentation time constants. In fact the loss of 1 mm 
grains is so great that the production factor for the whole 1 cm 
layer shown in figure 14 becomes more effected by the unexposed 
grains below the mixed zone and begins to decrease with time, 
as the mixed zone is depleted in 1 mm grains. 

When our presumed saturation time is also included (curves 

C) the production functions become well separated by grain size 

in a more surface-correlated manner. Note that the production 

function for 1 mm grains can reach an effective saturation due 

to fragmentation alone and all grain size fractions can reach 

saturation within typical residence times in the presence of 

our assumed saturation time T » 

8 


0.01 MT 
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For processes requiring only residence within the top 1 

governed by production functions and 

fragaentstion losses and grain saturation have such less effect 
except for soli particles whose dlaneter approaches or exceeds 
tne depth of the production zone, which Is 1 m in this case* 

We emphasize that saturation of a production function 
represents an average over all grains In the sawple and does 
not raean that all grains are saturated individually* In the 
less-nature soils modeled here at least half of the grains 
remain unexposed at the surface, others will reach saturation 
on the exposed portion of their surfaces. 

In figure 17 we have suzaaiarlzed the apparent grain size 
dependence D_ " implied by production functions G*(T_„)^* 

The .oUd lines give r. - 

the dashed lines give n « 1°8 jqIG*(Tj^) where 

X refers to EX and MM* Since all grains tend to reside equally 

within the first 1 mm the functions G* show strong surface 

rin 

correlation (n •» 1). The production functions for space 

exposure tend to be more volume correlated (n » 0) for short 

times. As time goes on both fragmentation and saturation cause 

G* to become more surface correlated, and this effect begins 
hX 

earlier for larger grains as suggested by Criswell (1975). 

These results apply to the early in-sltu development of less- 
mature soil. Mature soil will have been subjected to additional 
equilibrating effects of agglutination, comminution and mixing. 

To obtain the results described through figure 17 we have 
assumed a certain set of fragmentation and saturation time 
constants. However if tlu* time scale changes due to a variation 


^ figure 16), 
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In microoeteorold flux or If the eoler vind flux verles, then 
the effects of fragmentation and saturation as functions of time 
may be different* For example if for a given sample the appro- 
priate tlD» scale should be represented by 1 Mt ■ 10 tty due to a 
factor of 10 decrease in the sdcroaeteoroid flux, then would 
maintain the same value in units of MT (less f ragm«atatlon also) 
but would remain 0.01 My corresponding to 0*001 MT* For 
100-micron grains this would produce the curve D in figure IS, 
saturating at a factor of 10 below the corresponding curve C for 
1 MT ■> 1 My* Curves A and B would remain the same in this case* 
However if the tio« scale were to change to 1 MT ■ 10 My 
due to a reduction in crater size in a soil of greater cohesion, 
then both the fragmentation and saturation time constants will 
decrease by a factor of 10 in units of MT* This will produce 
the situation shown in figure 18 for the effective production 
function G* (T^)_ in the mixed zone of a 10 cm layer, where 
Tj^ ■ 0*03, 0.02, and 0.01 MT for 10-mlcron, lOO-micron, and 1 
mm grains, respectively, and T^ " 0*001 MT. In this case all 
grain size fractions can reach an effective saturation during 
typical residence times due to fragmentation alone* Note that 
for 1 MT - 10 My the typical residence times become 10 times 
longer than for 1 MT - 1 My* These properties may be useful in 
detecting time variations or variations in soil mechanical 
properties* In this context note also that ell track density 
parameters are proportional to the model time unit MT* 

In correlating track results with production functions 
we must keep in mind that the minimum, quartile, and median 
track densities are strongly controlled by the less irradiated 
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grains in a sample while Che production factors are usually 
controlled by the more irradiated grains. *lRierefore in a sample 
which mixes grains with the allowed wide range of local surface 
history there may be considerable variability in a correlation 
between tracks and noble gas concentration, for example. For a 
statistically average sample however the track parameters marked 
MM in figure 11 using the time axis should correspond to 
production factors given in figures 13 and 14 for 100~micron 
grains, while the curves marked CM in figure 11 using the T_„ 
time axis should corresponc' to the production functions in 
figures 15, 16, and 18 for 100-micron diameter grains. 

It is beyond the scope of this paper to compare our results 
in detail with the siany measured soil samples. Judging from the 
Shorty Crater data (Crosaz, 1979) and other track measurements 
the track parameters derived here seem reasonable for the postu- 
lated time scale 1 MT * 1 My. The situation for solar wind 
concentration is more complicated, being dependent on various 
mineral- and species-specific effects (Pilllnger, 1979). Our 
production factors seem reasonable in terms of average gas 
concentrations, however as Uieler e^ £l* (1980) have pointed out 
there is some difficulty when the distribution for individual 
grains is considered. Kirsten e^ (1971) have also measured 
gas concentrations in individual grains. 

Less-mature soils derived here have fewer exposed grains 
than expected for mature soils. However, a gas -poor soil 60051 
was found by Wieler e£ al . (1980) to contain some noble gas 
concentration In all 83 measured grains, apparently implying 
universal exposure for less-mature soils as well, contrary to 
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our model. These authors also measured track densities in 64 

grains from the same sample and found a median central density 
7 -2 

of less than 10 cm and only 19X with densities greater 
8 ~2 

than 10 cm . T)»se values are consistent with our model track 
results but inconsistent with universal space exposure* Wieler 
£t al . (1980) su^ert that grain fragmentation is respmislble 
for the lack of track density gradients which they also observet 
but fragments which lack such gradients have ne>wr been exposed 
and should lack SW gas concentration as well. 

It may therefore be necessary to find a mechanism yielding 
a surface concentration of SW elements in most grains without 
all being exposed to space* Such a mechanism could involve the 
Intermediate action of small grains* Criswell and Basu (1978) 
show that the size fraction less than 4 microns is generally 
capable of covering all larger grains which, together with the 
Intergrain forces acting on soil grains (Housley, 1980) will 
make it hard to avoid partial covering of most larger grains on 
the surface. In their study of microcraters on soil grains 
Poupeau et jd. (1973, 1978) find that even grains from mature 
soils have pits only on a small portion of their surfaces. 

If smaller grains coat part the surfaces of otherwise expo-- 
sed larger grains then they will acquire much of the SU gas. 
During subsequent burial through the uppermost few millimeters 
they will be subject repeatedly to shock compression by micro- 
meteoroid Impacts (StHff ler e^ , 1973) which dissipate a 
considerable fraction of their energy collapsing pore spaces 
between the soil grains (Braslau, 1970). Even weak shock can 
result in soil compactlSn and grain deformation, with minor 
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melting (Schaal and HUra, 1980). It seena plauaibla that during 
this process gas-loaded fine material coating soma of the larger 
grains could impart a concentration of SW atoma to the remaining 
unexpoaed grains residing near the surface. The process may be- 
come masked by subsequent exposure effects as the soil matures. 
The problem requires more work on the nature of grain-interface 
processes during soil c<»pactlon and involves a better under- 
standing of the depth profiles of SW elements in individual 
grains (£.£*• MtilJer et al«. 1976; Zinner et al.. 1978; Warhaut 
et el., 1979). 

If grains can acquire a SW concentration while residing 
within only the top millimeter of soil then it is clear from 
figure 17 that surface correlation can be obtained even for an 
Immature soil which is incompletely exposed. In addition if 
larger grains are subject to enlianced shielding due to inter- 
grain forces while residing on the surface then their exposure 
area controlled by the parameter f in equations 2 and 3 will be 
much less. Derived production functions G and G* will then be 

A A 

proportional to f < 1 and individual grains can reach microscop- 
ic saturation (Wieler ^ al. , 1980) at a total concentration 
much less than that for fully-exposed grains. Since the inter- 
grain forces primarily control smaller grains (i 10 microns, 
apparently) the smaller value of f should still be fairly 
constant for larger grains and increase again only for the 
smaller grains, as argued qualitatively by Housley (1980). 

In suonary, we have developed a model which allows us to 
derive distributions of track density, surface exposure times, 
and residence times within various zones as functions of grain 
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size, sample size and depth, and resid«ice tine. The results 
apply to samples which have resided within the top 10 cm for one 
time only and thus refer to relatively les8*^ture soils. Older 
surface layers which have been quickly burled deep enough to 
preserve their surface history undisturbed will also be describe 
ed by these results. Data from Shorty Crater has been used to 
help calibrate the model time scale, yielding reasonable track 
distributions and average exposure times. Important questions 
roiain concerning the distribution of solar wind elements among 
the grains in less-mature samples and concerning the role of 
very fine soil Mterial in controlling greln-slse dependences 
and solar wind gas distribution. This work can be extended to 
apply to mature soils which may have resided several times 
within the top 10 cm. 
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Figure Caption* 

1. A acheaatlc repraaentation of depth conaa daflnad by HESS* 2 
for a aodel regime. The model layer for a given regime cor* 
reapmidB to the surface aone for the next larger alee regiM* 

2. The ejecta origin i&ap developed for HESS.2* and 
are the central depth and rim radius* reepectively* of the 
excavated volume of the crater* h and r are the initial depth 
rnd radial distance, respectively* of a given soil particle* 
and r^ is the final radial distance of the ejected particle. 

3* A correlation between central track density and surface 

exposure time derived by M£Si^*2 for grains of diaiaeter 100 

microns after a typical residence time within the top 

1 esn of soil. The diagonal and dashed lines irilcate allowed 

lower and uppe*' limits, respectively* to the track density and 

the thin rectangle represents grains not directly exposed. 

In all figures 1 MT is a model time unit moat likely equal 
6 

to 10 years (1 My), unless otherwise noted. 

4. Typical cumulative distributions of space»exposure time t 
derived by MESS. 2 for grains of diameter 100 microns after a 
residence time within the top 1 mm of soil* for those 
grains which have been exposed to space (t > 0). 

5. Track parameters derived by MESS. 2 characterising tracks 
accumulated In the center of 100*mlcron diaiMter grains during 
residence within the top 1 mm of 5>uil* as functions of that 
residence time. MIN, Q* MED* 3Q* and MAX refer to minimum* 
first quartile, median* third quartlle* and maximum track 
densities* respectively. 
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6* Typical ct^ulative distributions of non-zsro rosidetics 
times within the top I nm of soil derived by ^^S•2 for 
grains of diameter 100 microns after a residence time T_., 
within the top 1 cm. 

7. Track parao^ters derived by MESS. 2 characterizing tracks 
accumulated In the center of 100*^cron diameter grains i^le 
residing at depths between 1 mm and 1 cm below the soil 
surface, as functions of grain residence time within the top 

1 cm. The dashed line gives the median track density 
accumulated at depths between 1 mm and 1 cm for those grains 
which have also resided within the top 1 mm» and hence repre- 
sents a partial density to be added to the track densities 
shown in figure 5. 

8. Typical cumulative distributions of non-zero residence 

times T within the top 1 cm of soil derived by HESS. 2 for 
CM 

grains of diameter 100 microns after a residence time T^„ 

I)M 

within the top 10 cm. 

9* Some properties of the mixed zone of a 10 cm surface layer 
derived by HESS. 2 after a typical layer residence time Tj^. 

The lower portion shows the cumulative distribution for depths 
of the mixed zone, collected over an area of radltis 25 cm. 

The upper portion is a scatter diagram for 100-micron diameter 
grains showing their depth plotted against their surface-zone 
residence time T^j (total tine spent within the top 1 cm). 

No general correlation of T^^^^ with depth is evident. 

10. The time development of the median depth of the mixed zone 
derived by MESS. 2 for certain regimes. The solid curves 
labelled MM, CM, and DM refer to 1 mm, 1 cm, and 10 cm thick 
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surface layers, respectively, which have locally survived 
destruction by relatively large impacts for residence times 

\ * ^MM* ^CIl* ^DM* dashed curve 

includes all sites, without the survival constraint* 

11. The minimum, first quartile, and Indian central track 
densities derived by MESS *2 for a sample of 100**ioicron dia*- 
meter grains as functions of residence time within the shal- 
lowest depth regime in which the entire sample has resided. 

Track con tribut long’s from greater depths are included according 
to derived residence time distributions. The time scales are 
lined up so that average values will correspond. 

12. Contours of equal quartile track density derived by MESS. 2 
for 100-micron diameter grains in a statistically average 
layer as functions of depth and layer residence time. The 
dashed curve gives the median depth of the mixed zone as in 
figure 10. A data point for the top of double drive tube 
74002 from Shorty Crater (circled dot) helps to calibrate 
the model time unit 1 MT as 10^ years (1 Ily). 

13. Effective area-time-per-mass production factors for surface 
exposure derived by MESS. 2 for typical grain diameters as 
functions of residence time within the top 1 cm* These curves 
are applicable to samples taken from the few-millimeters thick 
mixed zone of the 1 cm layer. Li figures 13 - 18 all curves 
labelled A do not include grain fragmentation or surface 
saturation, curves labelled B Include a loss by fragmentation 
due to direct impacts, curves labelled C include fragmentation 
plus an effective exposure tine corrected for a grain-surface 


saturation effect. 
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14. Effective area-tlne’-per'iaase production factora for surface 
exposure as in figure 13 except that these curves apply to a 
sample taken from the entire 1 cm surface layer. 

15* Effective area-tiine-*per-i&a88 production factors for surface 
exposure derived by MESS. 2 as functions of residence time 
within the top 10 cm. These curves apply to samples taken 
from the fe%/-centlmeters thick mixed zone of the 10 cm layer. 
The curve labelled D replaces curve C for D « 100 microns 
In the case of a factor of 10 decrease In the micrometeorold 
flux; curves A and B are unchanged in this case but the prob- 
able value of the model time unit 1 MT would become 10 tty., 

16. Effective area*'time-per'*mas8 production factors for 
residence within the top 1 mm derived by MESS. 2 as functions 
of residence tlTO within the top 10 cm. Ttese curves apply to 
samples taken from the mixed zone of the 10 cm layer. 

Effective residence times within 1 mm are much lees controlled 
by surface fragmentation and grain saturation except for soil 
particles whose diameter approaches or exceeds the depth of 
the assumed production zone. 

17. Assuming that the area-time-per--mass production factors 
follow a grain size dependence where n may be slowly 
varying, we approximate n In two size Intervals, 10 microns to 
100 microns (solid curves), and 100 microns to 1 mm (dashed 
curves) from the data shown in figures 15 and 16, as f.>nctlons 
of residence tin» within the top 10 cm. As before, A refers 
to no fragmentation or saturation, B Includes fragmentation 
loss at the surface, and C Includes frp 'dentation plus an 
effective exposure time due to grain saturation. 
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18. Effective area''tirae‘‘per-niaB8 pruduction factors for surface 
exposure as in figure 15 except that here the probable value 
of the model time unit 1 MT becomes 10 My due to a reduction 
in crater sise by a factor of about 2, in a soil with greater 
cohesion, for example. Both fragmentation and saturrtion 
will strongly effect all grain sizes in this situation* 
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